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FOREWORD 


This  is  Volume  III  of  three  separately  bound  volumes  in  which  are  reported  the 
research  completed  under  the  general  terms  of  the  Office  of  Civil  Defense  Subtask 
Number  4113E,  "Sensitivity  Analysis  of  Civil  Defense  Systems  and  Components." 

The  research  of  the  authors  was  very  ably  supported  by  Mr.  Herbert  Hill,  and 
Miss  Mary  B.  Woodside.  Mr.  Hill  assisted  in  the  development  of  the  models  and 
Miss  Woodside  performed  the  calculations  reported  in  several  of  the  Appendixes. 
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ABSTRACT 


This  sensitivity  analysis  employs  mathematical  models  which  estimate  the  total 
dose.  Maximum  ERD,  and  probability  of  casualty  or  fatality  for  an  individual  exposed 
to  a  particular  radiological  environment.  The  objective  of  the  analysis  is  to 
determine  the  relative  importance  of  the  parameters  normally  employed  in  evaluating 
Civil  Defense  systems;  i.e.,  to  determine  those  parameters  which,  through  large 
variance  or  inaccurate  estimates,  will  contribute  most  to  erroneous  evaluations  of 
CD  systems.  The  sensitivity  of  the  dose  (or  probability  of  casualty)  to  variations 
in  the  input  parameters  defining  the  radiological  environment  is  examined.  The 
total  Cone  model  is  analytical  and  the  ERD  model  is  computerized.  The  parameters 
examined  are  reference  intensity,  time  of  arrival,  time  outside  in  fallout,  radia¬ 
tion  decay  exponent,  ERD  recovery  fraction,  ERD  recovery  rate,  duration  of  fallout 
buildup,  and  protection  factor.  Sensitivity  indices  are  calculated  for  each  parame¬ 
ter.  The  sensitivity  index  is  defined  as  the  fractional  change  in  dose  (or  proba¬ 
bility  of  casualty  or  fatality)  divided  by  the  corresponding  fractional  change  in 
the  input  variable.  It  is  concluded  that  dose  and  casualty  computations  are  quite 
sensitive  to  errors  in  the  field  decay  exponent,  and  they  remain  sensitive  over 
the  examined  range.  Sensitivity  to  variations  in  fallout  reference  intensity 
and  protection  factor  are  high  over  the  whole  range  of  parameter  values.  Sensitivity 
of  time  of  arrival  of  fallout  can  be  quite  high  for  the  lower  values  of  the  parame¬ 
ter.  Sensitivity  of  dose  and  casualty  computations  to  the  remaining  parameters  is 
low  in  most  cases  of  interest.  Expansion  of  the  sensitivity  analysis  to  include 
parameters  other  than  fallout,  which  define  the  total  casualties  from  a  given  attack 
on  the  United  States,  is  necessary  before  further  conclusions  concerning  a  national 
vulnerability  analysis  can  be  drawn. 
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A  Generalized  Sensitivity  Analysis  of  CD  System 


I.  INTRODUCTION 

A.  Purpose 

An  analysis  of  the  effectiveness  of  the  civil  defense  shelter  system  must  take 
into  account  many  variables.  These  variables  include  some  that  determine  the  fall¬ 
out  environment  (fallout  reference  intensity,  its  time  of  arrival,  its  rate  of  build¬ 
up,  and  the  radioactivity  decay  rate);  some  that  determine  the  protection  against 
radioactivity  (the  gamma  ray  shielding  characteristics  of  available  shelter,  and  the 
time  required  for  travel  through  fallout  before  shelter  is  reached);  and  some  that 
influence  physiological  condition  (recovery  rate,  and  the  fraction  of  radiation 
damage  that  can  be  repaired).  The  research  reported  in  this  volume  is  directed  at 
evaluating  the  importance  of  uncertainty  or  inaccuracy  in  each  of  these  variables 
in  estimating  effectiveness  of  a  shelter  system  in  a  fallout  environment.  Effective¬ 
ness  is  measured  in  terms  of  Equivalent  Residual  Dose  (ERD) ,  probability  of  fatality, 
and  time  required  in  primary  shelter  for  the  ERD  to  return  to  a  specified  level. 

The  results  from  such  an  analysis  can  be  applied  in  two  ways:  (1)  determination 
of  research  priorities  and  requirements  for  increasing  the  validity  of  shelter 
systems  analysis,  and  (2)  improvement  of  the  models  or  formulations  that  are  used  for 
systems  analysis. 

B.  Procedure 

1.  General  -  - - - - - - - 

The  procedure  followed  in  this  volume  is  known  as  sensitivity  analysis, 
which  is  a  systematic  determination  of  the  effects  of  small  changes  in  input 
variables.  Such  an  analysis  requires  the  utilization  of  an  equation  or  a 
mathematical  model  that  relates  the  input  and  output  variables--such  as  reference 
intensity  and  estimates  of  casualties.  Each  input  variable  is  varied,  whilt 
all  others  are  kept  constant,  so  that  the  effects  of  these  variations  on  the 
output  are  measured  and  compared. 
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2. 


// 


Specific 


Because  of  the  number  and  nature  of  the  variables  that  influence  estimates 
of  the  effects  of  the  radioactive  fallout,  the  establishment  of  a  model  for 
sensitivity  analysis  is  complex.  The  major  steps  that  were  followed  in  this 
analysis  are  given  below  and  are  covered  in  more  detail  in  Section  II,  which 
follows.  A  generalized  flow  chart  of  the  procedure  is  shown  in  Figure  1, 
Research  Steps  in  Sensitivity  Analysis. 

Step  I.  Selection  of  Parameters  and  Values. 

Several  values  for  each  of  eight  factors  contributing  to  radiological 
effects  were  selected  as  input  parameters.  Combinations  of  these  values  were 
used  later  in  equations  (derived  while  performing  Step  II)  to  provide  the  cases 
upon  which  the  analysis  was  based.  (See  Section  II. A  below.) 

Step  II.  Derivation  of  Basic  Equations  to  Represent  Fallout  Effects. 


Equations  were  derived  for  determining  ERD  and  corresponding  probability 
of  fatality  (P^)  on  the  basis  of  the  eight  selected  factors.  (See  Section 
II. B  below.) 

Step  III.  Selection  of  Formulas  for  Determining  Relative  Sensitivity  of 
Input  Variables. 

The  following  relationship  was  established  as  the  sensitivity  index: 


Where  S^  is  the  sensitivity  index;  Y  is  the  measure  of  effectiveness  for  a 
selected  set  of  cases;  aY  is  an  arbitrary  value  change  in  Y;  X^,  is  the  value 
of  any  single  input  parameter  in  the  same  set  of  cases;  and  AX^  is  the  change 
in  this  input  necessary  to  cause  AY.  (See  Section  II. C  below.) 


Although  it  would  have  been  desirable  to  relate  the  equations  in  Step  II 
directly  to  the  sensitivity  ratio  in  Step  III,  it  was  not  possible  to  do  this 
readily  without  making  over-simplifying  assumptions.  Consequently,  intermediate 
steps  were  required. 


Step  IV.  Compilation  of  Cases  for  Analysis. 

Several  levels  of  each  input  variable  were  selected  and  the  Maximum  ERD 
and  probability  of  fatality  were  computed  for  each  combination  of  input 
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variables.  Using  a  CDC  1604  computer  and  the  Mainline  Program  which  was  designed 
specifically  for  this  purpose,  approximately  43,000  cases  were  developed.  After 
those  cases  were  excluded  in  which  the  Maximum  ERD  was  less  than  lOOr  or  greater 
than  700r,  a  total  of  10,143  cases  were  left  for  the  sensitivity  analysis. 

(See  Section  II. D  below.)  / 

Step  V.  Grouping  the  Cases  and  Calculation  of  Mean  Values. 

•  A  tape  containing  the  10,143  cases  was  transferred  to  a  general  purpose 
data  processing  program  written  for  the  IBM  7072  [Reference  1]  in  order  to  group 
the  selected  cases  into  subsets  for  analysis.  (See  Section  II. E  below.) 

Step  VI.  Determination  of  Model  Approximation  Equation. 

Using  the  10,143  selected  cases  and  multiple  regression,  the  same  general 
purpose  program  [Reference  1]  was  utilized  to  determine  a  single  equation 
(the  principal  model  approximation)  for  the  calculation  of  Maximum  ERD  as  a 
function  of  the  eight  input  variables.  (See  Section  II. F  below.) 

Step  VII.  Determination  of  Principal  Sensitivity  Indices. 

Ihe  effect  of  each  input  parameter  on  estimating  fatalities  was  calculated 
using  the  sensitivity  formula  (Equation  1)  stated  in  Step  III,  the  mean  values 
of  the  input  variables  determined  in  Step  V,  and  the  principal  model  approxima¬ 
tion  referred  to  in  Step  VI.  (See  Section  II. G  below.) 

Step  VIII.  Verification  of  the  Principal  Sensitivity  Indices. 

The  change  in  sensitivity  indices  was  examined  throughout  a  range  of  values 
using  a  Total  Dose  Model.  (See  Section  II. H  below.)  These  results  were  com¬ 
pared  to  those  contained  using  the  model  approximation  equation  of  Step  VI. 

(See  Section  il.I  below.) 

Step  IX.  Sensitivity  Analysis  of  the  Effect  of  Input  Parameters  on  Shelter 
Stay  Time. 

A  subgroup  of  cases  was  used  in  an  analysis  of  sensitivity  in  which  the 
dependent  variable  was  the  shelter  stay  time  required  for  the  ERD  to  return 
to  80r.  (See  Section  II. J  below.) 

Step  X.  Evaluation  of  Significance  of  Individual  Sensitivity  Indices. 

The  sensitivity  indices  of  the  various  parameters  were  then  evaluated  in 
view  of  the  accuracy  with  which  these  parameters  could  be  expected  to  be 
established  as  inputs  to  any  analysis  of  system  effectiveness.  (See  Section 
II. K  below.) 
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II.  THE  SENSITIVITY  ANALYSIS 


A.  Selection  of  Parameters  and  Values 

In  order  to  obtain  a  representative  sample  of  cases  on  which  to  base  the  analysis 
of  sensitivity,  a  range  of  reasonable  values  was  assigned  to  each  of  the  input  variable 
commonly  used  in  shelter  systems  analysis.  These  values  are  summarized  in  Table  I.  It 
will  be  noted  that  in  all  cases  except  that  for  protection  factor  (PF) ,  the  values  are 
evenly  spaced. 


TABLE  I 

Values  and  Symbols  for  the  Input  Parameters 
Used  in  the  Sensitivity  Analysis 


Parameter 

Symbol 

Values 

Reference  Intensity  (r/hr) 

I 

o 

300,  1500,  2700,  3900 

Time  of  Arrival  (hrs) 

ta 

1,  4,  7 

Buildup  Duration  Factor 

E 

1.13,  3.64,  6.15 

Protection  Factor 

PF 

2,  10,  40,  100,  500 

Time  in  Fallout  (hrs) 

T2 

0.0,  0.3,  0.6 

Field  Decay  Exponent 

Z 

-1.0,  -1.2,  -1.4 

ERD  Recovery  Fraction 

F 

0.85,  0.90,  0.95 

ERD  Recovery  Rate  (Fraction/day) 

B 

0.020,  0.025,  0.030 

Fallout  hazard  (Risk)  data  on  200  Standard  Metropolitan  Statistical  Areas 
(SMSA's)  were  used  as  the  basis  for  selection  of  values  of  reference  intensity  and 
time  of  arrival.  These  data  were  extracted  from  an  OCD  source.  Times  of  arrival 
between  one  hour  and  seven  hours  occurred  in  about  50  percent  of  those  cases; 
reference  intensities  were  between  300  r/hr  and  3900  r/hr  for  94  percent  of  the  cases. 
(See  Volume  II,  Figs.  31  and  32.) 

The  buildup  duration  factor  (E)  is  a  factor  by  which  time  of  ariival  is  multiplied 
in  order  to  represent  estimates  of  possible  time  between  the  beginning  and  cessation 
of  fallout.  Appendix  A  contains  a  discussion  of  buildup  functions. 

Ihe  protection  factor  (PF)  values  were  based  partially  on  National  Fallout  Shelter 
Survey  data  summarized  in  Reference  [3].  Identified  shelter  spaces  with  PF's  of  40 
and  above  could  accommodate  roughly  20  percent  of  the  United  States  population.  PF's 
of  40,  100,  and  500  were  selected  to  represent  the  three  levels  of  PF's  between  40 
and  99,  between  100  and  250,  and  greater  than  250,  respectively.  Factors  of  10  md  2 
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were  used  as  representatives  of  the  PF  of  population  in  residential  basements  and  in 
homes  without  basements,  respectively.  ; 

/ime  in  fallout,  prior  to  shelter  occupancy  was  set  arbitrarily  at  0.0,  0.3, 
and  0.6  hours. 

In  determining  equivalent  residual  dose  (ERD) ,  values  used  for  the  radiation 
field  decay  exponent  and  for  the  ERD  recovery  fraction  and  ERD  recovery  rate  were 
varied  about  the  generally  accepted  values  used  in  systems  analysis. 


B.  Model  Basic  Equations 

The  following  basic  equations  constitute  the  analytical  formulation  of  the 
model.  It  was  upon  these  equations  that  the  Mainline  Computer  Program  (see  Appendix 
B)  was  developed. 

1.  Equivalent  Residual  Dose  (ERD) 


ERD  (n)  ■  D  +  d 
'  n  n 


Dn  “  (1"B)  °n-l  +  F  Cl  *(t)  dt 

dn  *  dn-l  +  G-F)  Cl  *(C)  dt 


(1) 

(2) 

(3) 


where : 


D  ■  Recoverable  dose  remaining  during  time  interval  n-1  to  n 
n 


dose  remaining  during  time  interval  n-1  to  n 

I 

n1 

e  of  recovery  for  reparable  damage  (fraction/ 


(roentgens). 

d  ■  Nonrecoverabld 
n 

(roentgens) . 

B  ■  Fractional  rat] 

— — - unit  time). 

F  *  Fraction  of  dose  which  is  reparable. 

(1-F)  *  Fraction  of  dose  which  is  nonreparable. 

n  ■  Number  of  time  units  after  detonation,  n  *  1,  2,  3,  ...» 
when  increments  (n)  in  time  are  24-hour  days.  (Time  is 

incremented  in  minutes  for  the  first  24  hours  and  in  days 

thereafter.) 

t  *  Number  of  hours  after  detonation. 

$(t)  ■  Field  Intensity  in  r/hr  at  time  (t)  after  detonation. 


and : 


t(t)  -  0 


for:  t  <  T. 

**  A 


where : 


sensitivity  index  for  each  input  parameter. 


We  begin  the  mathematical  development  by  first  defining  the  sensitivity  index 
as  the  ratio  of  the  percent  change  in  the  output,  or  dependent  variable,  to  the  per 
cent  change  in  an  input  parameter  [Reference  4], 

Thus,  the  sensitivity  index  (S^)  of  the  iC^  input  parameter,  measured  from 
mean  values,  is: 

[a.  gc  in  the  dependent  variable  (AY)] 

c  [Mean  .  ’ue  in  the  dependent  variable  (Y)] 

5i  - - 

[Change  in  the  1  input  parameter  (AXj,)] 

[Mean  value  of  the  iC^  input  parameter  (X^)] 


or: 


(4) 


D.  Calculation  of  ERD  and  Probability  of  Fatality 

The  model  basic  equations  were  used  with  each  combination  of  input  values  in 
the  Mainline  Program^  on  the  CDC  1604  computer. 


For  each  combination  of  input  parameters,  the  time  required  for  the  buildup 
of  the  radiation  field  is  calculated.  Next,  the  dose  a  person  would  receive,  if 
any,  while  outside  the  shelter  and  the  dose  at  the  end  of  the  first  day  inside  is 
computed  and  totaled.  The  dose  is  then  split  into  a  recoverable  e^d  nonrecoverable 
fraction  and  a  test  is  made  to  determine  if  the  probability  of  fatality  equals 
"one."  If  so,  the  program  returns  to  the  initial  state  and  selects  another  input 
combi  .tion. 

The  ERD  is  accumulated  and  updated  by  recoverable  and  nonrecoverable  fractions 
and  the  following  tests  are  made  on  the  daily  increments: 

(1)  If  a  Maximum  ERD  has  been  reached,  its  values  and  date  of  occurrence  are 
recorded  for  output. 


—  An  additional  independent  variable,  D, ,  was  included  to  set  three  possible  levels 
to  which  the  ERD  could  return — 80r,  50r,  and  20r.  There  was  a  total  of  42,740 
cases  for  all  combinations  of  input  values. 
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(2)  If  a  Maximum  ERD  has  been  reached  and  the  nonrecoverable  ERD  has  not  ex¬ 
ceeded  the  allowable  dose,  and  the  accumulated  ERD  is  less  than  the 
allowable  ERD,  then  the  day  is  recorded  for  output. 

Eventually  the  test  conditions  will  be  met  and  the  program  will  branch  out  of 
the  day-by-day  calculation  of  the  ERD.  Once  out  of  this  routine,  the  program  will 
take  the  Maximum  ERD  value  and  calculate  the  probability  of  fatality  and  other  de¬ 
pendent  variables.  _  _ ___ _ _ _ 

Hie  dependent  variables  calculated  by  this  program  are: 

(1)  Maximum  ERD. 

(2)  Probability  of  fatality. 

(3)  Probability  of  being  a  noneffective. 

(4)  Day  in  which  the  Maximum  ERD  occurs. 

(5)  Day  in  which  allowable  ERD  occurs. 

(6)  Nonrecoverable  dose  on  that  day. 

Cases  were  excluded  in  which  the  Maximum  ERD  was  grearer  than  700r  cr  less  than 
lOOr,  thereby  leaving  for  analysis  a  group  of  10,143  cases  (Gi.oup  0  of  Table  II, 
Basic  and  Subgroup  Divisions  of  the  Mainline  Program  Cases).  Thus,  the  sensitivity 
analysis  was  based  on  cases  in  which  the  interaction  of  input  parameter  values  re¬ 
sulted  in  situations  of  maximum  interest--that  is,  cases  in  which  even  small  varia¬ 
tions  in  input  parameters  could  affect  both  casualties  and  fatalities. 

E.  Grouping  of  Cases  and  Calculations  of  Mean  Values 

Subgroups  of  Group  0  (which  contained  all  cases  between  lOOr  anc  700r  Maximum 
ERD)  were  established  for  each  of  the  four  input  values  of  reference  intensity.  In 
addition,  it  was  desirable  to  combine  two  of  these  subgroups  in  order  to  have  a 
subgroup  with  "moderately"  high  reference  intensity  (Subgroup  20  of  Table  II). 

In  order  to  examine  the  significance  of  the  time-in-shelter  data,  another 
basic  group  of  data  was  established  (Group  1).  Into  this  group  were  put  cases  in 
which  the  ERD  would  eventually  reach  80r  or  less,  because  80r  is  considered  to  be 
the  threshold  below  which  physical  disability  is  unlikely  to  exist.  Within  this 
latter  basic  group,  subgroups  based  on  reference  intensity,  allowable  ERD,  and 
protection  factor  were  established.  Table  II  shows  the  basic  and  subgroup  divisions. 

The  cases  were  transferred  from  the  Mainline  Program  by  a  tape  to  the  IBM  7072 
computer.  A  general  purpose  data  processing  program  (TSAR)  at  Duke  University  was 
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TABLE  II 


Baste  and  Subgroup  Divisions  of  the  Mainline  Program  Cases 


BASIC  CROUPS  SUB' 

■CROUPS 

NUMBER 

SORTING  PARAMETER 

OF  CASES 

Group  0 

10,143 

Reference  Intensity  (r/hr) 

(parameter  Combinations 
with  : 

5 

2,115 

I  -  300 

lOOr  <  Max.  ERD  <  700r) 

6 

3,108 

1°  -  1500 

7 

3,431 

1°  -  2700 

8 

1,48? 

1°  -  3900 

20 

of 

(sum 

7  &  8) 

4,920 

1°  -  3063 
o 

Group  1 

3,481 

(Cases  which  return  to 
an  allowable  ERD  <  80r) 

Reference  Intensity  (r/hr) 

9 

909 

I  -  300 

10 

1,041 

1°  »  1500 

11 

1,058 

1°  -  2700 

12 

473 

1°  -  3900 

0 

Allowable  ERD  (r) 

2 

2,069 

D  «  80 

3 

1,204 

d£  -  50 

4 

208 

d£  -  20 

Protection 

13 

987 

PF  *  2 

14 

1,291 

PF  -  10 

15 

795 

PF  *  40 

16 

321 

PF  «  100 

■ 

17 

87 

PF  -  500 

Multiple  Sorting  Parameters 

18 

352 

D4  -  80r 

0  <  P0  <  .5 

2  — 

t2  -  0 

19 

240 

D.  -  50^ 

4 

0  <  P„  <  .5 

2  - 

t2  -  0 

-  10 

used  for  the  grouping  of  data  into  sets  rRe£ 


The  same  program  was  used  for  determining  the  mean  values  of  the  input  parameters 


for  the  basic  groups  of  cases  and  for  each  of] 
of  the  Basic  Parameters,  shows  the  pertinent 
groups. 


F.  Determination  of  the  Principal  Model  Approximation 


rence  1 


the  subgroups.  Table  III,  Mean  Values 
heans  for  each  of  these  groups  and  sub- 


For  purposes  of  computation  it  was  desirable  to  have  a  single  equation  that 

del.  Using  a  stepwise  regression  to 
,  logarithmic,  and  exponential)  were 
fitted  to  the  sets  of  data  produced  by  that  nlodel.  This  was  accomplished  also  by 


approximates  the  Mainline  Program  computer  nx> 
obtain  this,  three  basic  functions  (quadratic 


the  TSAR  general  purpose  program.  Each  funct 


«[ Reference  11.  |, Aquation  5  best  fitted  the  dajta  and  was  made  the  principal  model 
approximation 


ion  was  fitted  by  th 2  step  procedure 


f  (X)  =>  Max.  ERD  ■ 

e(-2. 7468)  j  (.932  +  .0001499PF 


-  .0446F  +  5.354B  -  .006493E) 


_  (.4346  -  . 7954Z)  p_(-1.0949  +  /04181T,  4-  .06268Z) 
TA  “  i  2 

(-3.074  -  . 1161E)  (-1.1553) 


(5) 


G. 


I 

The  statistical  characteristics  of  the  principal  model  approximation  equation 


are: 


(1)  Multiple  correlation  coefficient: 


0.9991. 


(2)  Standard  error  of  the  estimate:  0.02421  (i.e.,  c  of  ERD  +  6r). 


(3)  Value  of  F-ratio:  886.785  (a  further  indication  of  the  high  significance 
of  the  equation). 

(4)  Data  Group  0  (Table  II). 


Calculation  of  the  Sensitivity  Indices 

We  now  have  a  formula  (Equation  5)  relating  the  input  variables  to  the  dependent 
variable  Maximum  ERD,  and  a  formula  for  the  sensitivity  index  (Equation  4).  To 
calculate  the  sensitivity  indices  at  the  mearjs.  Equation  4  is  placed  in  differential 
form: 


— ^  Identification  of  the  variables  is  in  Table  I. 
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Individuals  in  500  PF  shelters  only  reached  or  exceeded  the  Maximum  ERD  if  they  receive  appreciable  dose  outside. 
Their  time-outside  average  was  twice  any  other  group. 


Let: 

so; 


-  X  (Y)"1 

&xi  Xi  W 

V  -  f(X) 


S1 


(*> 


(6) 


Note  that  f(X)  is  the  principal  model  approximation  (See  Equatioti  5): 

^1  ^2  ^3  *"4  ^5 

f  (X)  •  K  I  T.  PF  Z  B 

O  A 


(7) 


where : 

.  C1  .932  +  .0001499PF  -  .0446F  +  5.354B  -  .006493E 

C2  -  .4346  -  . 7954Z 

C3  -  -1.0949  +  .04181T2  +  .06268Z 

C4  -  -3.074  -  . 1161E 

C5  *  -1.1553 
-2.7468 

K  *  e 


Ihe  partial  derivatives  of  Equation  5  with  respect  to  each  independent  variable 
(Xt)  are: 


Ci"1  C,  C-  C,  C, 

.  k  c  i  1  T.  2  PF  3  Z  4  B  5 

c»T  1  o  A 

o 

mffXV  Ci  C2“l  C-  C4  C5 

'  K  C2  TA  "  2  B 

zim.  .  K  K  PF03’1  ZC‘  b'5 

3PF  K2  o  *A  "  Z 

TA  PF  2  B 

im  .  K  K.  .1  Cl  T.°2  Pf'3  ZV‘  B°5 

az  4  o  A 


him 

bF 


>KKj 
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-  14  - 


F  -  .9 
B  -  .025 
E  -  2.57 


and  using  the  Equation  4,  we  calculate  the  Sensitivity  Indices  (S^). 

TABLE  IV 

Principal  Sensitivity  Indices 


Index 

Variable  Denoted  by  Index  Subscript 

Value 

sz 

Field  Decay  Exponent 

-4.22 

SPF 

Protection  Factor 

-0.96 

SI 

o 

Reference  Intensity 

+1.02 

Sta 

Time  of  Arrival 

-0.50 

SF 

ERD  Recovery  Fraction 

-0.30 

SB 

ERD  Recovery  Rate 

-0.14 

SE 

Buildup  Duration  Factor 

-0.18 

ST 

l2 

Time  in  Fallout 

40.05 

It  is  upon  these  values  that  the  conclusions  and  recommendations  will  be 
developed. 


The  above  sensitivity  indices  are  obtained  for  Maximum  ERD  as  the  dependent 
variable.  However,  over  the  dose  range  where  probability  of  casualty  is  neither 
0  nor  1,  che  relation  between  probability  sensitivity  (PSp  and  dose  sensitivity 
(S^)  Is  given  by  the  relation: 

PSi  "  Si  (1  *  p>  <»> 

where  p  ■  probability  of  casualty  «  a  •  Dose  +  b.  Hence,  the  relative  rank  of  the 
Sf  remains  the  same  when  evaluated  at  a  given  set  of  input  variables. 

H.  Total  Dose  Model 

The  major  computational  difficulties  in  determining  the  above  sensitivity 
indices  arise  because  of  the  use  of  equivalent  residual  dose  instead  of  total  dose 
as  a  measure  of  biological  effect.  This  section  will  describe  a  simpler  model  which 


-  15 


ft 

/ 

/ 

allows  the  calculation  of  the  sensitivity  indices  analytically.  In  this  model  the 
dose  rate  buildup  function  is  assumed  to  be  linear,  and  total  dose  is  used  as  the 
dependent  parameter. 


A  new  parameter,  To,  is  introduced  in  this  approach.  Tq  is  the  time  period  from 


entrance  into  shelter  until  the  cutoff  of  the  total  dose  calculations.  Ihus, 

(T^  +  T2  ■]+  Tq)  represents  the  complete  time  period  from  time  of  detonation  until  com¬ 


pletion  of  the  total  dose  calculation.  Using  the  parameters  previously  defined  (see 


/ 


Table  I)  and  the  new  parameter  T^,  the  total  dose  D  is  given  by  the  expression: 


n  Xo 
D‘P? 


2TaZ+1  (l+E)~E 


( 


<PP-1>  T22  + 


(1 


(Ta  +  T2  +  T0>1’2  -  Ta1_Z  (1+E) 


"1 


The  sensitivity  indices  can  then  be  calculated  directly 

(PF-1)  T2 


s  T2  Xo  T2 

s__  *  ■!.  1  •  T  ■  — ■  ■ 


•12  -  SI2 


D  PF 


XAZ+1  (1+E)‘E 


z  +  <TA  +  X2  +  V 


ta  sta  d 


DJPF  T 


Z-l 


TaZ  (Ta  +  T2  +  To)'Z  -  (1+E) 1_z 


2 (1+E) 


i_  [(pp-i)e«)(|)  ♦  o-i)  e2 


c  PF  m  o 

rF  "  SPF  *  D  2DPF 


2  „2„  2. 


TA  <1+E>  e 


(T,  -  E‘T,) 


(1-Z)  [ 


TA  +  T2  +  To)1_Z  -  TaX“Z  (1+E) 


i 


ST  *  If  '  k  *  (TA  +  T2  +  V 

o  o 


-2 


(10) 


(ii) 


(12) 


(13) 


(14) 
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(15) 


_  M3  E 
SE  "  ?E  '  D 


I  E 
o 


D  PF  (1+E) 


Z+l 


(1-PFVH-E  di-2)  2 

7J.1  1  1  7 


2T  z+1  E2 
A 


+  ^4^  ta1-z  -  TaW  (1+E) 


(16) 


m>  z  I02 

SZ  "  'Z  ’  D  “  D  PF 


E 


2  -i 

V  ET  7 

(pf-1)  2tT  +  “2^  ln  ta(1+e>  •ta(i+e);‘  [-1] 

A  J 


4 


(Ta  +  T2  +  To)1‘Z  -  [Ta(1+F)1 


‘1 


(1-Z) 

(ife  |<Ta  +  T2  +  To>‘'Z  >"  <TA  +  T2  +  V 


-  ;ta(i+e);1_z  in  [ta(i+e) 


(17) 


The  total  dose  sensitivity  indices  (S ^ )  can  be  readily  converted  to  probability 
of  casualty  sensitivity  indices  (PS ^) »  if  desired.  In  the  region  of  interest 
(probability  neither  zero  nor  one),  the  dose  is  related  to  probability  of  casualty 
(p)  by  a  linear  relation: 


p  *  aD  +  b. 


(18) 


Hence : 


X. 

i 

P 


si  *  jb 


(19) 


Only  dose  sensitivity  will  be  discussed  in  the  remainder  of  this  analysis,  since  the 
relative  rank  o'  both  S  and  PS  are  the  same  for  a  given  set  of  values  of  the  inde¬ 
pendent  variables. 


In  Table  V,  Comparison  of  Sensitivity  Indices,  the  values  of  the  S^  obtained 
from  Equations  7  through  13  by  evaluation  at  the  mean  values  of  the  variables  listed 
in  section  G  above  are  compared  to  the  values  of  S^  listed  in  Table  IV  (calculated 
from  the  ERD  model). 
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Index  Variable  Denoted  by  Subscript  ERD  Model  Value 


otal  Dose  Model  Value 


Field  Decay  Exponent 


Protection  Factor 


Reference  Intensity 
Time  of  Arrival 

! 

Buildup  Duration  Factor 


-4.22 

-0.96 


+1.02 


-0.50 


-4.36 


-0.99 

+1.00 

-0.50 

-0.25 


Time  in  Fallout 


+0.05 


+0.02 


In  Figures  2  thru  6,  the  values  of  the  sensitivity  indices  are  plotted  as  a 
function  of  the  indexed  variable.  While  the.  rest  of  the  independent  variables  are 
held  constant  at  the  mean  values  given  previously,  two  values  for  PF  (39.1  and  500) 
and  two  values  of  Tq  (96  and  240  hours)  are  plotted  for  each  curve.  Figure  7  shows 
the  increase  in  the  sensitivity  index  for  (time  in  fallout)  as  Z  (the  decay  rate) 
increases.  Figure  8  shows  the  dose  sensitivity  to  T  as  a  function  of  T  for  the 

OiO 

two  protection  factors.  This  curve  indicates  that  total  dose  is  not  particularly 

'  0  •  * 

sensitive  to  Tq  when  Tq  is  measured  after  100-200  hours,  especially  for  high  PF 
shelters.  In  fact,  where  Maximum  ERD  occurs  within  the  first  jfew  days,  it  can  be 
expected  that  the  results  from  the  total  dose  sensitivity  analysis  model  will  be 

essentially  equivalent  to  the  results  from  the  ERD  model. 

V  1  . 

I.  Validity  of  Ranking  the  Sensitivity  Indices 


In  Tables  IV  and  V,  the  sensitivity  indices  were  evaluated  for  one  set  of  input 

i 

variables.  However,  the  relative  values  of  these  Sj,  might  change  when  evaluated  for 
radiological  environments  other  than  the  one  defined  by  the  mean  values  in  section  G. 
To  investigate  this  problem,  the  following  conclusions  were  drawn  based  on  Figures 
2  through  8  and  other  calculations  of  sensitivity  indices: 


(1)  Sj  ,  sensitivity  index  for  reference  intensity  (I0)>  jts  1  over  the  entire 
range  of  values;  thus,  Iq  remains  a  relatively  important  parameter  in  all 


(2)  Sn_  is  essentially  constant  at  1  for  most  cases  of  interest;  thus,  protec- 
PF 

tion  factor  (PF)  is  a  relatively  important  parameter  in  all  cases.  (See 
Figure  2) . 
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Protection  Factor  (pp) 

Fig.  2.  Sensitivity  Index  of  PF  as  a  Function  of  PF 


Time  of  Arrival  (T.)  in  Hours 
A 


12  3  4  5  6  7 


Fig.  4.  Sensitivity  Index  (S_)  as  a  Function  of  Buildup  Factor  (E) 


Sensitivity  Index  (S„,  )  as  a  Function  of  Time  in  Fallout 


Fig.  6.  Sensitivity  Index  (S„)  as  a  Function  of  Decay  Constant  (Z) 


SensitivityJLndex  for  Time  in  Fallout  (S  )  as  a  Function  of  Decay  Constant 


me  (hrs)  over  which  Total  Dose  is  Calculated 


(3)  S  remains  high  for  all  cases,  thus  Z  is  extremely  critical  in  vulnerability 

M 

analyses.  (See  Figure  6). 

(A)  For  analysis  observing  only  Maximum  ERD  reached  while  in  primary  shelter, 

ERD  concepts  and  a  total  dose  approximation  of  Maximum  ERD  yield  essentially 
the  same  results  in  sensitivity  analysis.  (This  does  not  imply  that  they 
yield  the  same  results  in  casualty  calculations.) 

(5)  can  become  quite  large  for  short  times  of  arrival  of  fallout  (T^)  and 
high  protection  factors.  In  areas  where  this  combination  of  environments 

exists,  T^  must  be  treated  carefully  in  vulnerability  analyses.  (See  . — 

Figure  3). 

(6)  In  most  cases  of  interest,  the  other  sensitivity  indices  remain  appreciably 
less  than  1,  though  some  change  significantly  as  the  parameter  values  are 
varied. 

Analysis  of  Time  Required  in  Shelter 
1.  Procedure 

In  the  preceding  analysis,  the  independent  variable  was  Maximum  ERD.  It 
is  also  important  to  know  how  the  input  variables  may  affect  the  time  required 
in  shelter  after  the  beginning  of  fallout. 

In  the  analysis  of  time  required  in  shelter,  only  those  cases  from  the 
Mainline  Program  were  considered  in  which  the  probability  of  a  noneffective 
(P2)  was  greater  than  zero  and  less  than  0.5.  Thus,  the  largest  Maximum  ERD 
to  be  considered  was  200r  and  the  smallest  Maximum  ERD  was  lOOr.  (See  the  dose 
response  relationships  in  Section  II. B. 2.)  Another  condition  was  that  the  ERD 
returned  to  80r.  It  was  assumed  that  the  reduction  of  ERD  to  this  level  would 
mark  approximately  the  end  of  the  continuous  stay  in  shelters.  Time  in  fallout 
outside  shelter  was  also  excluded  from  consideration. 

The  TSAR  general  purpose  program  was  used  to  set  up  subgroups  18  and  19 
(see  Table  II)  and  to  develop  an  equation  by  multiple  regression.  The  equation 
given  below  determines  the  number  of  days  in  shelter  until  ERD  returns  to  80r. 

When:  lOOr  £  Maximum  ERD  <  200r, 


(20) 


f(X)*  -  N,  **  4.47  -  44Z  +  2.26  +  1.6D,  -  1.55(Z  x  -£) 

&  r  *  J  r  r 

+  926 (Z  x  3)  *  .7920 (D^  x  F)  -  19.8(D3  x  B) 

* 

Multiple  Correlation  Coefficient:  0.9717 
Value  of  F-ratio:  830 


Where : 

N2  ■  Duration  of  shelter  phase  until  ERD  falls  to  80r 
Z  ■  Field  decay  exponent 

Jo  ■  Reference  intensity/protection  factor 
PF 

«*  Maximum  ERD  prior  to 

B  ■  ERD  recovery  rate 
F  «  ERD  recovery  fraction 

The  mean  value  of  stay  time  for  the  cases  considered  (Y  in  the  sensitivity 
formula)  was  found  to  be  47  days. 

The  sensitivity  indices  were  calculated  using  Equation  20,  the  differentia, 
form  of  Equation  4: 


CL 


e  ,  moo  i 

i  ?Xi  ’  47 

witn : 

*  Y  *  f(X)  (see  Equation  20), 

N2  *  Y  =  47  days. 


(21) 


2.  Results 

The  sensitivity  indices  for  the  input  parameters  are  shown  in  Table  VI, 
Sensitivity  Analysis  of  Time  Required  in  Shelter.  It  is  seen  that  changes  in 
the  decay  exponent  (Z)  are  the  most  significant  in  affecting  the  length  of  stay 
in  shelter. 
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TABLE  VI 


// 

/ 
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Sensitivity  Analysis  of  Time  Required  In  Shelter 


Parameter 

1 

Mean  Value 

*i 

Gradient 

Sensitivity  Index 
(Equation  21) 

Si 

Z 

1.260 

-  144.000 

-3.860 

F 

0.914 

-  144.000 

-2.220 

°3 

144.000 

0.373 

1.140 

B 

0.025 

-1700.000 

-0.918 

RP 

79.600 

0.331 

0.560 

The  independent  variables  (which  affect  the  length  of  stay  in  shelter)  are 
ranked  on  the  sensitivity  index  in  decreasing  importance  as  follows: 


1.  Field  decay  exponent,  Z  ■  3.9 

2.  ERD  recovery  fraction,  F  »  2.2 

3.  Maximum  ERD,  *  1.1 

4.  ERD  recovery  rate,  B  »  0.9 


5.  Reference  intensity/protection  factor  ratio,  RP  *•  0.6 

Of  equal  interest  to  the  parameter  sensitivity  analysis  is  the  actual 

timing  of  Maximum  ERD  and  the  number  of  days'  stay  in  the  primary  shelter  before 

recovery  to  an  allowable  dose  is  accomplished.  For  a  shelter  PF  of  73,  the 

Maximum  ERD  of  144r  occurred  on  the  7th  day  after  the  attack.  On  the  average, 

3/ 

47  days  elapsed  before  the  ERD  recovered  to  the  allowable  dose  level  of  80r.— 
Figure  9  shows  the  frequency  of  occurrence  plotted  against  the  duration  in  days 
before  reaching  an  ERD  of  80r.  It  can  be  seen  that  there  are  cases  in  which 
stay  time  would  exceed  100  days  if  the  ERD  of  shelter  occupants  is  to  return 
to  80r . 

K.  Evaluation  of  Significance  of  Sensitivity  Indices  in  Shelter  Systems  Analysis 
1.  General 

The  objective  of  shelter  systems  analysis  is  to  predict  as  accurately  as 
possible  the  adequacy  of  existing  or  potential  protective  measures  against 

3/ 

—  If  instead,  an  allowable  level  of  50r  (Subgroup  19)  were  chosen,  then  the  average 
length  of  stay  would  increase  by  33  days  (from  47  to  80). 
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Days  for  Return  of  ERD  to  80r  (N2) 

Fig.  9.  Frequency  of  Occurrence  of  Stay  Times  vs  Days  for  Return  of  ERD 
(for  Cases  in  Which  Max  ERD  <  200) 


fallout.  In  general,  the  better  the  parameters  that  determine  the  adequacy  are 
known,  the  better  the  predictions  will  be. 

The  sensitivity  indices  derived  in  the  above  analysis  give  an  indication  of 
the  effects  of  not  knowing  the  various  input  parameters  exactly;  the  analysis 
does  not  yet  take  into  account  how  accurately  these  parameters  might  be  known. 

In  the  following  sections,  an  attempt  will  be  made  to  identify  the  approxi¬ 
mate  magnitude  of  uncertainty  to  be  expected  in  the  eight  parameters  considered. 

2.  Reference  Intensity  and  Time  of  Ar-ival 

Of  the  eight  parameters  that  were  considered  in  the  sensitivity  analysis, 
reference  intensity  and  time  of  arrival  are  largely  dependent  upon  the  choice 
of  strategy  in  the  attack;  consequently,  these  parameters  cannot  be  known 
accurately.  Desired  ground  zeros  and  the  occurrence  of  air  or  surface  burst  are 
the  choice  of  the  attacker. 

At  any  given  point,  the  reference  intensity  and  time  of  arrival  are  also 
affected  by  wind  speed  and  direction. 

The  RISK  type  programs  are  attempts  to  represent  the  combined  effects  of 
strategic  decisions  and  wind  conditions.  In  the  choice  of  input  variables  for 
the  sensitivity  analysis  (section  II. A),  it  was  seen  that  reference  intensities 
between  300  and  3900r/hr  and  times  of  arrival  between  1  and  7  hours  seem  to 
represent  the  range  of  intensities  with  which  fallout-only  shelter  systems  are 
primarily  concerned.  Variations  from  the  means  of  these  input  parameters  of 
+  75%  would  seem  to  be  a  representative  estimate  of  the  probeble  range  of  values 
for  reference  intensity  and  time  of  arrival  in  an  analysis  of  a  shelter  system. 

3.  Buildup  Factor 

The  values  of  the  buildup  multiplying  factor  used  in  the  sensitivity 
analysis  ranged  from  1.13  to  6.15.  These  extreme  values  vary  from  their  mean 
by  about  707.  and  this  is  proposed  as  a  range  of  values  that  could  be  expected 
in  systems  analysis. 

4.  Protection  Factors 

The  National  Fallout  Shelter  Survey  (NFSS)  produced  data  on  the  present 
number  and  the  protection  factors  of  shelter  spaces.  However,  there  are  various 
reasons  why  these  data  alone  may  not  represent  the  actual  shelter  posture  of  the 
population. 
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First,  based  on  a  limited  sample  of  shelter  [Reference  5],  it  appears  that 
procedures  used  in  the  Phase  1  shelter  survey  generally  underestimated  the 
actual  protection  factor  by  a  considerable  amount.  The  current  updating  of 
NFSS  data  is  gradually  improving  the  accuracy  of  information  on  shelter  status, 
but  a  complete  reevaluation  of  all  NFSS  shelters  is  not  contemplated  by  OCD. 

At  the  present  time  the  average  potential  protection  factor  of  the  United  States 
population  may  be  considerably  more  than  presently  represented  in  Phase  1  or 
Phase  2  data.  However,  ever,  if  the  protect  tur.  factors  w\.re  accurately  deter¬ 
mined,  there  still  would  remain  uncertainties  as  to  shelter  utilization  that 
would  influence  the  results  of  shelter  analysis.  The  question  of  shelter 
utilization  was  presented  in  Volume  II  where  it  was  considered  in  terms  of 
actual  SMSA's.  (See  Volume  II,  A  Sensitivity  Analysis  of  Selected  Parameters 
Based  on  8  SMSA's.l 

Consequently,  random  variations  of  807.  about  the  mean  can  be  considered 
reasonable  in  estimating  the  significance  of  the  sensitivity  index  for  protec¬ 
tion  factor. 

5.  Time  in  Fallout 

Estimates  of  length  of  time  outside  in  fallout  are  arbitrary.  By  assuming 
adequate  warning,  no  one  would  need  to  be  outside.  However,  there  may  also  be 
cases  where  large  numbers  of  the  population  are  in  fallout  for  periods  ex¬ 
ceeding  the  0.6  hours  arbitrarily  used  in  the  sensitivity  analysis. 


Although  the  sensitivity  index  for  time  in  fallout  is  small,  the  importance 
of  this  parameter  in  systems  analysis  must  be  determined  in  conjunction  with 
actual  locations  and  categories  of  shelters.  The  combined  effects  of  these  are 
analyzed  in  Volume  II. 

6.  Decay  Rate  . - . — - - - - — .  . - 

Ihe  decay  rate  may  vary  with  the  construction  of  the  weapons  and  to  a  lesser 
degree  with  the  nature  of  the  surface  materials  at  ground  zero.  Variations 
probably  will  not  be  more  than  plus  or  minus  25  percent  of  the  1.2  rate  of 
decay,  which  is  based  on  empirical  data. 

7.  Biological  Factors 

The  fraction  of  the  dose  that  is  reparable  (F)  and  the  rate  of  recovery 
(B)  are  important  parts  of  the  Blair  Formula  for  calculating  ERD.  However 
according  to  Reference  [6],  a  firm  experimental  basis  appears  to  be  lacking 
for  the  values  presently  used.  The  author  in  the  same  reference  also  questions 


whether  the  constants  are  independent  of  time. 

In  this  study,  the  average  value  of  F,  the  recovery  fraction,  is  taken  as 
0.90;  variations  of  10  percent  would  seem  likely.  Ihe  average  value  of  B,  the 
recovery  rate,  is  0.025  and  probably  could  vary  by  70  percent. 

8.  Comparison 

We  then  may  compare  the  sensitivity  indices  with  the  rough  estimates  of  the 
accuracy  to  be  expected  in  input  variables.  These  are  sunmarized  in  Table  VII. 
An  approximation  of  the  relative  importance  of  each  parameter  in  contributing  to 
total  uncertainty  in  shelter  analysis  is  obtained  by  multiplying  the  range  of 
uncertainty  by  the  sensitivity  indices. 

TABLE  VII 

t- 

Comparison  of  Sensitivity  Indices  and  Uncertainty  in  Input  Parameters 


Variable 

S  (ERD 
Model) 

S  (Total 

Dose  Model) 

Uncertainty 

Approximate 

Relative 

Importance 

Field  Decay  Exponent 

-4.22 

-4.36 

25% 

1.09 

Refererce  Intensity 

+1.02 

+1.00 

75% 

0.75 

Protection  Factor 

-0.96 

-0.99 

80% 

0.80 

Time  of  Arrival 

-0.50 

-0.50 

75% 

0.35 

Buildup  Duration  Factor 

-0.18 

-0.25 

70% 

0.18 

ERD  Recovery  Rate 

-0.14 

--- 

70% 

0.10 

Time  Outside  in  Fallout 

+0.05 

+0.02 

100% 

0.05 

ERD  Recovery  Fraction 

-0.30 

--- 

10% 

0.03 

CONCLUSIONS  AND  RECOMMENDATIONS 


A.  Conclusions 


Considering  both  the  sensitivity  of  each  parameter  and  the  range  of  uncertainty 
that  each  is  likely  to  have,  it  is  possible  to  draw  some  conclusions  as  to  their 
effect  on  the  results  of  systems  analysis. 

In  analyses  for  estimates  of  Maximum  ERD,  the  following  have  the  most  signi¬ 
ficance:  Field  decay  exponent  (Z),  reference  intensity  (Iq),  and  protection  factor 
(PF).  Having  less  significance  are  time  of  arrival  (T^) ,  ERD  recovery  fraction  (F), 
ERD  recovery  rate  (B),  and  buildup  duration  factor  (E). 


Time  in  fallout  (T2)  has  a  low  sensitivity  index.  However,  extending  the  range 
of  values  that  were  assumed  for  it  beyond  0.6  hours  may  have  a  considerable  effect 
on  the  magnitude  of  this  number.  The  importance  of  time  in  fallout  was  better  shown 
in  Volume  II,  where  it  was  analyzed  as  a  factor  in  shelter  utilization. 

In  analysis  of  required  time  in  shelter,  field  decay  exponent  (2)  and  ERD  re¬ 
covery  rate  (B)  are  the  most  significant. 

Of  the  above  parameters,  reference  intensity  and  time  of  arrival  are  almost 
entirely  dependent  upon  the  strategic  choice  of  the  attacker  and  wind  conditions. 
Consequently,  better  estimates  of  their  effects  in  systems  analysis  is  dependent 
upon  military  intelligence  more  than  upon  additional  research  and  analysis. 

The  importance  of  the  field  decay  exponent  in  determining  both  fatalities  and 
shelter  stay  makes  it  important  to  have  continuing  analysis  of  the  validity  of  the 
1.2  value.  There  may  be  appreciable  gains  in  accuracy  possible  since  the  exponent 
is  primarily  dependent  upon  the  physics  of  the  radioactive  decay  process. 

Increased  knowledge  of  the  shelter  protection  factors  available  to  the  population 
is  important  for  the  making  of  valid  systems  analyses.  The  effect  of  underestimates 
of  protection  factors  can  be  taken  into  account  by  using  a  multiplying  factor. 

However,  it  is  important  to  define  the  factor  more  closely  than  is  possible  now. 

This  would  not  necessarily  require  large  expendif.res  for  obtaining  new  data  or  for 
correcting  data  that  now  exist.  Statistical  analysis  based  on  the  parts  of  the 
survey  that  have  been  updated  or  examined  in  detail  could  provide  a  much  better 
estimate  than  evidently  now  exists  of  the  overall  shelter  that  would  be  available  to 
the  population. 

Identification  of  the  biological  recovery  process  is  particularly  important  in 
estimating  time  required  in  shelter.  Possible  approaches  to  reducing  tne  range  of 
uncertainty  of  this  input  are  given  in  Reference  [6], 

Finally,  efforts  to  define  more  closely  the  limits  of  any  input  variables  should 
be  viewed  in  the  light  of  their  relative  effect  on  output  data.  For  example,  attempts 
to  define  protection  factors  accurately  may  make  a  change  in  the  relative  importance 
of  that  parameter  that  is  small  compared  to  the  relative  importance  of  reference 
intensity,  which  cannot  be  accurately  predicted  by  any  research  means. 

The  sensitivity  index  of  the  buildup  duration  factor  (E)  is  based  on  "moderate" 
conditions  where  the  time  of  arrival  approximates  3.3  hours  and  (E)  approximates  2.6. 

In  observing  the  value  of  the  sensitivity  index  S^,,  it  was  found  that  the  index  is 
quite  dependent  on  changes  in  (E).  Some  further  analysis  of  the  effect  of  the  build¬ 
up  duration  factor  is  discussed  in  Appendix  A. 
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The  analysis  based  upon  return  of  ERD  to  80r  (the  threshold  of  appearance  of 
clinically  detectable  symptoms  of  radiation  sickness),  showed  that  shelter  stay  times 
of  much  more  than  two  weeks  could  be  necessary.  It  should  be  noted  that  this  analysis 
included  only  cases  in  which  the  Maximum  ERD  was  200r  or  less.  Without  these  restric¬ 
tions,  even  longer  required  stay  times  would  result.  This  is  very  significant  in 
planning  for  the  stocking  of  shelters,  recovery  planning,  and  for  the  establishment 
of  policies  for  postattack  management  of  supplies. 

In  summary  the  conclusions  from  the  sensitivity  analysis  are: 

(1)  Casualty  calculations  are  quite  sensitive  to  errors  in  the  field  decay 
exponent.  The  sensitivity  of  dose  or  casualties  to  the  field  decay  rate 
remains  at  approximately  -4.0  over  the  range  examined.  (The  decay  rate 
was  varied  from  Z  «*  1.1  to  Z  •*  1.4.) 

(2)  Sensitivity  to  variations  in  fallout  reference  intensity  and  protection 
factor  are  high. over  the  whole  range  of  parameter  values.  Hence,  precise 
knowledge  of  the  fallout  shelter  posture  and  the  fallout  reference  inten¬ 
sity  is  much  more  essential  to  accurate  vulnerability  analysis  than  the 
remaining  parameters  in  most  cases  of  interest. 

(3)  Sensitivity  to  time  of  arrival  of  fallout  can  be  quite  high,  in  some 
radiological  environments. 

(4)  Sensitivity  of  casualty  computations  to  the  remaining  parameters  is  low 
in  most  cases  of  interest. 

(5)  Expansion  of  the  sensitivity  analysis  to  include  parameters  other  than 
fallout,  which  define  the  total  casualties  from  a  given  attack  on  the 
United  States,  is  necessary  before  further  conclusions  concerning  a 
national  vulnerability  analysis  can  be  drawn. 

B.  Recommendations 

1.  Vulnerability  analyses  should  employ  protection  factors  computed  by  best 
available  methods,  and  research  and/or  surveys  to  improve  protection  factor  data 
should  be  encouraged. 

2.  Because  oi  the  sensitivity  of  systems  analysis  results  to  the  field  decay 

-1  2 

exponent,  continuing  analysis  of  the  validity  of  the  t  ’  decay  law  should  be  made. 

3.  Additional  study  is  required  to  establish  the  sensitivity  of  fatalities, 
casualties,  and  dose  to  duration  of  shelter  stay. 
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4.  Sensitivity  analysis  should  be  extended  to  include  the  parameters  defining 
the  effects  of  blast  and  alternative  measures  of  effectiveness  (casualties  by  type, 
lose  distribution  of  survivors,  etc.). 


5.  Sensitivity  analysis  should  be  applied  to  identification  of  the  important 
tost/effectiveness  parameters  used  in  the  budget  allocation  model  (see  Volume  I, 
k  Cost/Effectiveness  Computer  Procedure  for  Optimum  Allocation  of  Fallout  Shelter 


System  Funds  Under  Uniform  or  Variable  Risk  Assumptions' 
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Appendix  A 


The  Intensity  Buildup  of  a  Radioactive  Fallout  Field 


In  calculating  the  radiatior  dose  received  at  a  point  in  the  fallout  field,  it 
is  necessary  to  determine  the  rate  at  which  the  fallout  particles  arc  deposited  from 
the1 cloud.  The  easiest  method  o i  representing  this  buildup  is  to  plot  the  fraction 
of  jtotal  fallout  mass  deposited  £,s  a  function  of  time.  If  is  the  time  the  first 


parades  arrive,  and  is  the  lime  the  last  particles  are  deposited,  then  the 

duration  of  the  buildup  period,  ’’  -  T  ,  is  approximated  by  a  relation  of  the  form 

v/  A 

tc'ta  =  ta‘e  (a_1> 

where  E  is  the  buildup  duration  factor  referred  to  in  this  volume. 

i  _ 

Figure  A-l  presents  the  buildup  function  ''Reference  A-l.  which  was  obtained 

from  weapons  tests.  In  addition^  two  approximations  to  this  function  are  plotted 

i 

in  ’Figure  A-l  for  =  1  and  E  ■; 2.5.  The  analytic  form  for  these  approximate 
curves  is  given  by  Equations  A-2  and  A-3. 


!<o  ■  v  ■('■=) 

-  i  -  cos  J-J* 


(A-2) 


(A-3) 


Equation  A-3  gives  a  somewhat  better  fit  of  the  curve  from  Reference  rA-l', 

’  |  * 
hence  this  buildup  function  was  used  in  the  computer  model  that  calculated  Maximum 

i 

ERD.  However,  the  tctal  dose  expressions  cannot  be  integrated  in  closed  form  if 
Equation  (A-3)  is  used,  hence  Eqbation  (A-2)  was  used  as  the  buildup  function  for  the 

analytical  total  dose  model.  As  might  be  expected  from  Figure  A-l,  little  difference 

' 

was  found  between  sensitivity  indices  from,  the  two  models. 
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Time  after  Detonation 


Fallout  Mass  Buildup  Functions 


Appendix  B 

_ The  Mainline  Model  and  Computer  Program 

I.  THE  OPERATIONAL  TRANSFORM  OR  MAINLINE  MODEL 


Attack  Environment  Formulation 


l.  Introduction 

Ihe  first  step  in  formulating  the  Mainline  Program  is  proper  characteriza¬ 
tion  of  the  attack  environment.  The  interrelation  between  the  initial  intensity, 
field  buildup,  the  field  decay,  and  the  field  attenuation  by  shelter,  are  the 
four  functions  that  predict  the  dose  rate.  To  determine  their  relative  im¬ 
portance,  it  becomes  mandatory  that  tnis  interrelation  is  accurately  formulated. 

The  underlying  assumptions  throughout  this  section  are: 


(a)  The  Reference  Intensity  of  the  radioactive  field  (in  roentgens  per 
hour)  will  conform  to  the  measurements  given  in  the  NREC  RISK-type 
data  (i.e.,  an  intensity  (Iq)  at  H  +  1  hours)  as  well  as  the  distribu¬ 
tion  of  intensity  and  arrival  of  fallout  obtained  also  from  NREC  RISK- 
type  data.  | 


(b)  Ihe  measure  of  protection  obtained  by  the  emploj'ment  of  shelters  will 
be  the  Protection  Factor  (PF)  obtained  from  "Thjj 
Phases  1  and  2. 


! 

Fallout  Shelter  Surveys 


2.  Decay  of  Field  Intensity 


Since  a  radioactive  particle  decays  with  time,  our  problem  is  to  determine 
the  appropriate  mathematical  representation,  when  applied  to  residual  fields 
used  in  Civil  Defense  Systems  Analysis.  The  widely  used  function  shown  in 
Figure  B-l  and  Equation  (B-l)  will  be  used. 


where 


0(t)  ■  Field  Decay  Function,  in  r/hr. 
IQ  «  Reference  Intensity,  in  r/hr. 
t  ■  Time  in  hours. 

Z  •  Empirical  Constant. 


Fig.  B-l.  Decay  of  Radioactive  Field 

The  two  important  constants  introduced  with  this  function  are:  the 
empirically  derived  exponent  (Z),  and  the  reference  intensity  (Iq).  The  refer¬ 
ence  intensity  is  one  of  the  two  interface  parameters  connecting  the  fallout 
shelter  analysis  model  with  the  '*war  games"  which  establish  the  probable  attack 
environment. 

3.  Buildup  of  Field  Intensity 
a.  Introduction 

The  problem  now  is  to  characterize  analytically  the  field  intensity 
during  the  buildup  phase.  The  fact  that  there  is  no  existing  data  within 
our  files  on  which  the  usual  methods  of  curve  fitting  can  be  applied,  forced 
us  to  use  a  more  general  form  of  analysis.  Our  analysis  uses  what  little 
significant  information  there  is  and  employs  as  much  intuition  as  seems 
necessary,  useful,  or  feasible.  Although  there  are  two  relevant  functions 


(B-2) 


and 


(B-3) 
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they  are  incompatible  over  a  wide  range  of  times  of  arrival.  The  relation 
time  of  peak  intensity  equals  twice  the  time  of  arrival--is  taken  as  being 
more  nearly  valid.  In  the  field  buildup,  two  distinct  phenomena  are 
present.  The  first  is  the  arrival  (at  ground  level)  of  the  radioactive 
particles.  The  second  is  the  independent  phenomena  of  the  radioactive 
decay  of  these  particles.  Thus,  characterizations  of  these  two  phenomena 
will  be  independent,  and  the  field  intensity  represented  by  the  product. 

The  hypothesis  implies  a  particle  buildup  function  of  the  following 
general  form  given  in  Figure  B-2,  which,  when  integrated,  yields  the  time 
dependent  arrival  function  given  in  Figure  B-3. 


Time  t  - ► 

Fig.  B-2.  Particle  Arrival  Rate 


Fig.  B-3.  Accumulation  of  Particles 
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The  produce  of  the  function*  represented  in  Figure  B-l  and  Figure  B-3 
would  give  the  field  intensity  shown  in  Figure  B-4. 


Fig.  B-4.  Field  Intensity  vs.  Time 

b.  Formulation 

From  several  general  forms  a  displaced  Cosine  function  was  chosen 

for  simplicity  and  compatibility  with  the  function  (t  ■  2T.)  to  represent 
^  P  A 

the  fallout  particles  accumulated  (Reference  Equation  B-4). 

F(t)  -  j  £  -  cos  ">]  (S’4) 

where 

T,  <  t  <  T 

A  -  -  C 

and 

1  ■  TC  •  V 

-z 

The  field  decay  function  (i.e.,  t  )  times  the  buildup  function  times  the 

reference  intensity  yields  the  final  representation  of  the  field  intensity 

over  time.  For  time  between  T.  and  T  : 

A  0 

!(t)  -  Ic  ~  [1  -  cos(i^-^  n)J  (B-5) 
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Thus,  the  characterization  of  the  fi<  Id  intensity  is  in  two  parts;  the 
first  during  accunulation--fron  T  (time  of  arrival)  to  Tc  (time  of 
cessation).  The  second,  after  all  fallout  has  occurred  (Reference  Equation 
B-l),  describes  the  field. 

In  order  to  comply  with  T  =  2T,,  the  value  of  must  be  fixed.  The 

p  A 

following  are  values  of  for  three  values  of  Z: 

Z  -  1.0  tj  •  1.13  3Ta 

Z  -  1.2  tv  -  1.175T. 


Z  -  1.4 


1.20T, 


Attenuation  of  Field  Intensity 

a.  The  Problem 

The  population  or  an  individual  may  move  through  a  series  of  shelter 
conditions  over  time.  The  duration  of  each  shelter  phase,  although 
nominally  under  operational  control,  has  well-defined  bounds  based  upon 
the  attack  environment.  The  initiation  of  the  recovery  phase  depends  upon 
minimizing  the  stay  in  shelter.  Thus,  the  influence  of  shelters  when 
modeled  must  take  into  account  the  operational  requirement  to  minimize 
shelter  stay  time. 

b.  Formulation 

Fallout  shelters  are  classified  by  a  protection  factor  (PF) .  The 
protection  factor  function  operates  on  the  unprotected  radiation  field 
to  yield  the  attenuated  field  that  exists  within  that  shelter. 

There  are  four  major  operational  phases  (Reference  Figure  B-5). 

These  are: 

(1)  An  initial  unsheltered  condition  which  corresponds  with  a  person 
moving  through  a  fallout  field  before  entering  shelter.  This 
phase  will  be  called  "time  outside."—^  It  is  related  to  the 
efficiency  of  warning  systems,  drills  and  training,  and  shelter 
assignment  plans. 

(2)  Phase  2  is  the  primary  shelter  phase.  In  most  current  models  it 
is  the  only  one  considered.  Operational  plans  usually  specify 


— ■  The  symbol  used  in  this  report  for  initial  time  outside  in  a  fallout  field  is  T£. 
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Field  Intensity  (r/hr) 


(3) 


Phase  3  is  operationally  assoc 
begin  recovery.  The  objective 


a  two-week  stay,  and  they  are  itocked  on  that  basis.  The  control¬ 
lable  variable  is  the  length  of  stay  in  the  phase.  Our  model 
considers  the  length  of  stay  dependent  upon  an  objective  of  the 
third  phase  and  is  subject  to  explicit  calculation. 


ated  with  emerging  from  shelter  to 
in  this  period  is  to  obtain  produc¬ 


tive  effort  (usually  at  a  PF  or  1)  without  any  physiological 
radiation  symptoms.  To  do  this  introduces  a  decision  parameter 
(an  allowable  ERD),  an  operational  parameter,  and  the  equivalent 
PF  required.  The  allowable  ERD  lies  between  the  limits  of 
clinical  effects  (25r)  and  symptomatic  effects  (lOOr). 


(4)  Tne  fourth  phase  is  the  return 


to  normal.  It  begins  when  the 


individual's  recovery  rate  equils  the  dose  received  in  the  un¬ 
sheltered  field.  It  is  characterized  by  the  lack  of  control 
variables,  decisions  required,  etc.  In  other  words,  the  effects 
of  normal  biological  recovery  processes  exceed  the  effects  of  the 
radiation  fields. 

The  problem  domain  investigated  in  this  volume  extends  only  to  Phase  3, 
time  (T^).  (Practical  considerations  precluded  extensive  investigation  of 
Phases  3  and  4  under  the  present  contract.) 


Fig.  B-5.  Attenuation  of  Field  Intensity 
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In  Fig.  B-5,  the  solid  line  represents  the  sheltered  point  intensity. 
In  Fig.  B-5,  the  dotted  line  represents  the  field  intensity  if: 


where 


PF^  «  PF2  *  PFj  "  1. 


Phase 

1 

PF1(t)  -  1,  0  <  t  <  T2 

Phase 

2 

PF2(t)  -  PF,  T2  <  t  <  T3 

Phase 

3 

PF3(t)  *  Equivalent  PF(t),  T3  <  t  <  T^ 

Phase 

4 

PF4(t)  -  1,  T4  <  t. 

5.  Summary  Attack  Environment  Model  ' .  .  . . . . — . —  - . -  — 

The  attack  environment  radioactive  field  intensity  at  any  time  (t)  is  the 
product  of:  the  Reference  Intensity  (IQ) ;  the  particle  buildup  function  (F(t)); 
the  field  decay  function  (•' (t)) ;  and  the  protection  factor  function  (PF(t)). 


*(t)  -  R  •  F(t) 


9(t) 

PF(t) 


(B-6) 


Or  by  substitution: 


when  T,  <  t  <  T,  +  -,T 

A  ""  A  A 


*(t)  -  0  when  t  <  TA  (B-7) 

1  - 

5(t)  -  ~  t  when  t  >  TA  +  rjT^ 


where: 

t(t)  *  Field  Intensity  in  r/hr  at  time  t  after  the  originating 
explosion  event. 

Iq  *  The  reference  intensity  in  r/hr,  one  hour  after  the  explosion 
event. 

t  ■  Time  in  hours  measured  from  H. 

T,  ■  Time  of  arrival  measured  from  H  hour. 

A 

PF,  ■  Protection  factor  at  time  t. 
k 

Z  ■  empirically  derived  field  decay  constant. 

n  ■  T„  -  T,  ■  fallout  interval. 

1  C  A 
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Biological  Response 


1.  Introduction 

Ihe  state,  efficiency,  and  outcome  of  a  shelter  system  analysis  is  measured 
by  the  biological  impact,  the  usual  criteria  being  the  number  of  fatalities  or 
lives  saved.  Fatalities  are  calculated  by  using  a  probability  of  response  to 
the  total  radiation  dose  received. 

At  the  present  state  of  understanding,  the  equivalent  residual  dose  model 
is  the  best  available  formulation  of  biological  response  to  dose  received. 

2.  Equivalent  Residual  Pose 

To  represent  analytically  the  absorbed  gamma  radiation  field  in  terms  of 
a  time  dependent  dose  "state"  requires  taking  into  account  the  observed  phenomena 
of  biological  recovery.  Ihe  ERD  model  relates  the  field  intensity  and  biological 
recovery  rate  to  give  a  single  "dose"  value,  approximating  the  physiological 
state  of  a  human. 

Ihe  ERD  at  any  point  in  time  consists  of  two  parts.  Ihe  nonrecover able 
portion  is  the  accumulative  amount  of  the  nonrecoverable  dose  received  over 
time.  Ihe  recoverable  portion  has  two  additive  parts;  one  is  a  fraction  of  the 
field  contribution  within  the  unit  time  considered,  and  the  second  is  the 
recoverable  dose  carried  over  from  the  previous  period.  The  basic  statement 

of  ERD  (recoverable  dose  (D  )  in  time  increment,  n)  is  the  sum  of  that  not  yet 

“  th 

received  from  previous  periods  and  the  new  contribution  during  the  n  time 

period. 

Dn  -  (l-B)Dn_1  +  F  /n  «(t)dt  (B-8) 

1  n-1 

where  • _ _ ; _ _ _ 

B  is  the  percent  per  unit  dose  recovered  in  period  n,  and  F  is  the 
recoverable  fraction. 

Nonrecoverable  dose  (d  )  in  time  increment  (n)  is  the  accumulative  sum  of  all 

n 

nonreparable  doses: 

dn  "  dn-l  +  (1'F)  /"  ®(0«»t.  (B-9) 

n-1 

Ihe  practical  problem  of  describing  the  function  of  the  ERD(n)  in  a  closed, 
well-behaved  form  requires  that  it  be  approximated  by  finite  differences. 

Finite  difference  problems,  which  require  recursive  solutions,  lead  to  computer 
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Implementation.  A  solution  is  to  use  the  basic  difference  equations  with  n 

in  one  day  intervals  and  with  d  being  computed  daily,  then  added  for  the 

^  th 

equivalent  residual  dose  on  the  nfc  day. 
where 

D  ■  Recoverable  dose  in  time  interval  n-1  to  n, 
n  ’ 

d  ■  Nonrecoverable  dose  in  time  interval  n-1  to  n, 
n  * 

I  «  Reference  intensity, 

F  “  Fraction  of  dose  which  does  reparable  damage, 

(1-F)  ■  Fraction  of  dose  which  does  nonreparable  damage, 

B  -  Daily  fraction  of  dose  repaired, 

■  Hie  time  of  arrival  of  the  first  measurable  radioactive  fallout, 

Z  *  Decay  constant, 

PF  *  Protection  factor. 

3*  Probability 

The  problem  is  now  to  relate  the  equivalent  residual  dose  which  an  individ 
ual  accumulates  to  the  probable  biological  state.  The  usual  specifications  for 
radiological  induced  states  are  fatalities  and  casulaties,  with  fatality  being 
a  type  of  casualty.  Because  this  causes  some  confusion,  we  will  introduce  the 
term  noneffective.  For  a  given  Maximum  ERD,  we  will  then  have  three  probabili¬ 
ties: 

(1)  Death. 

(2)  Noneffect ive--The  condition  of  being  alive  but  physically  unable  to 
perform  a  productive  task. 

(3)  Normal  and  Marginal  Effective — Individuals  with  acceptably  low  Overt 
physiological  symptoms  (below  lOOr). 

The  following  linear  functions  are  used  to  approximate  the  empirically  derived 
probability  of  casualty  vs.  ERD  curves  (reference  Figure  B-6). 

a.  Probability  of  Fatality 

Pj  -  0  when  ERD  <  200r 

Pj  _  ■QfeXj.~Sgil.200  when  200 r  <  ERD  <  700r 

Pj  -  1.  when  700r  <  ERD. 
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S  lOOr  200r  300r  400r  500r  600r  700r 

0* 


Equivalent  Residual  Dose  (ERD) 
Fig.  B-6.  Dose  Response 


Probability  of  Noneffectiveness 


when  ERL  <  lOOr 
when  lOOr  <  ERD  <  200r 
when  200r  <  ERD  <  300r 
when  300r  <  ERD  <  700r 
when  700r  <  ERD. 


(B-U) 


(B-12) 


C.  Mainline  Program 
1.  Introduction 

The  models,  previously  described,  are  combined  into  a  single  nK>del--a 
computer  program.  The  finite  difference  method,  particularly  adaptable  to 
computer  progTamning,  is  necessary  because  of  several  functional  discontinuities 
which  preclude  integrating  the  individual  models  simultaneously.  In  translating 
algebraic  formulas  into  FORTRAN  coding,  a  change  in  notation  is  required.  To 
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assist  in  identifying  corresponding  it  rw ,  a  Glossary  of  Syr.brls  use-1  in  the 
programming,  and  the  corresponding  symbols  found  in  mathematical  formulation, 
are  listed  in  Table  B-I. 


TABLE  B-I 
Program  Symbols 


Parameters 

Model  Symbols 

•Fr.ograrr  ing  Symbols 

Reference  Intensity 

V 

R 

Time  of  Arrival 

.  TA 

I/. 

Time  Outside 

T2 

. 

Allowable  Dose 

D4 

!>• 

Field  Decay  Constant 

Z 

Z 

Recovery  Fraction 

F 

Recovery  Rate 

ft 

0 

Buildup  Duration 

f:;a 

Protection  Factor 

PF 

PF 

Protection  Factor  Categories 

K 

Maximum  No.  PF  Categories 

KMAX 

Time  (the  variable  of 
integration) 

Hours 

Days 

V 

l 

n 

T 

Total  (ERD)  Dose 

ERD  (n) 

D 

Recoverable  D.-se 

D 

n 

D1 

Konrecover<iMe  Dose 

ct 

n 

D2 

Day  Maximm  EJID  Occurred 

%T 

Ml 

Field  Intercity 

:  ft) 

RDN  or  FLF 

Day  Allowable  Dose  Occurred 

N2 

N’2 

Probability  of  Fatality 

Fl 

PI 

Probability  of  Noneffective 

P2 

P2 

In  the  input  sections  of  the  program,  if  a  programing  symbol  given  in 
Table  3-1  is  followed  by  L,  H,  or  S,  the  nuirbvr  represented  by  the  symbol  is 
either  the  lowest,  highest,  01  Incremental  vala  taken  by  the  variable  repre¬ 
sented  by  the  preceding  letters  (i,t.,  TAL  means  the  lowest  value  the  time  of 
arrival  will  rake),  ,  , 
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71*  irt  ke/ed  t-v  '  (hm  diagram  by  the  organ  rat  lor  ti 

•  tat  raw  t  lire  russet  with  the  *:.'•■'■>  *:«*-*>*  .  Specifically,  each  of  the 
coding  lire  numbers  art  four  dig^  "  irst  digit  represents  the  major  routine 
being  calculated  («ee  Table  B-II);  t  second  and  third  digits  correspond  to  the 
flew  diagram  box  number  (xx) ;  the  fourth  or  units-position  digit  is  "open”  to 
allow  for  sequencing  corrections,  additions,  and  expansions  (0). 

TABLE  B-II 

_  ..  ...  _  Major  Program  Routines  _ _ _ _ _ _ _ _ _ _ _J„ 

Code 

Number  Series 


Input 

lxxO 

-  -  Iterations  &  Loops - — 

- - 2xx0 

Calculation  of  Buildup 

Duration  in  hours 

3xx0 

Initial  Conditions 

4xx0 

Calculate  1st  Day  Doses 

5xx0 

Day  by  Day  Doses 

6xx0 

Outcome  Probabilities 

7xx0 

Formats  &  ,fHousekeeping" 

Steps 

8xx0 

Explanation  of  the  Flow  Diagrams 

The  first  section  of  the  program  (block  numbers  1-36)  initializes  the 
input  data  and  sets  up  the  control  loops  for  the  variables.  Block  numbers  36 
through  40  select  the  appropriate  duration  multipliers  to  calculate  the  dura¬ 
tion  of  Particle  Buildup  (T3), 

T3  -  (1.133  +  ETA)*TA  when  2  =1. 

- - - — _  (1.175  +  ETA)*TA - when  2  -  1.2 - .  (B-13) 

■  (1.200  +  ETA)*TA  when  2  *  1.4 

the  range  of  Buildup  Duration  Factor  (ETA)  is  between  zero  (corresponding  to 
the  displaced  sine  function)  and  five  (corresponding  •  ' th  5TA  ^  -  TA  function). 

The  program  sets  the  initial  conditions  in  block  47.  Clock  Time  (T)  begins 
at  time  of  arrival  (TA)  and  the  Day  Counter  (N)  begins  at  one,  the  time  is 
incremented  by  3  minutes  or  0.05  hrs. 
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Field  ir.tcr'iity  (r  DV )  ter  unit  of  tine  ii  calculated  ir.  Mocks  4i,  *5, 
utinR  the  pr oduc t  of  Reference  Irtensitv  (r/hr).  Particle  Arrival  functicr,  a -6 
Decay  function. 

Limit  on  T3  is: 

RDN  -  |  |l  -  cos<2^)*JT*Z  0  <  <■  -  (B-14) 

The  test,  block  42,  determines  if  there  is  any  movement  outside  in  fallout--or 
if  T2  is  positive. 


Block  52  is  similar  to  47  except  the  field  is  attenuated  by  the  protection 
factor.  Block  53  tests  to  determine  if  buildup  of  particles  is  continuing. 
Block  54  increments  time  out  of  buildup  phase  so  that  field  intensity  is 
governed  only  by  the  decay  function  corresponding  to  the  algebraic  formula. 

0  -  T“Z  .  (B-15) 


Block  56  tests  to  determine  if  the  1st  day  is  completed.  Block  55  increments 
time  by  12  minutes  during  buildup.  Blocks  49,  50,  51  are  identical  to  44,  45, 
46  but  are  not  linked  because  of  a  change  in  the  time  intervals. 


FLD  in  block  59  is  equivalent  to  RDN  but  is  used  to  account  for  dose  when 
the  time  step  is  by  days. 


Ht) 


f2«+TA  t-Z  dt 


1*Z  |  25+TA 

C(t)  -  1  ^  ■  |T  when  Z  4  1  (B-16) 

e(t) . 

v  '  Z-l 


Block  60  does  not  have  a  time  term  since  the  integration  in  block  59  was  over 
a  finite  time. 


Blocks  64,  65,  66,  67  and  68  sets  up  transition  from  the  first  day,  to 
day-by-day  calculations.  At  blocks  o4-C3  "u''  Mm  ERD  division  into  recoverabl 
and  nonrecoverable  fractions  is  made.  Block  66  steps  time  by  12  hours.  Block  67 
sets  maximum  dose  register.  Block  69  solves  equation  for  Field  Intensity  summed 
over  24  hours  but  corrected  according  to  PF, 


__  R*T*Z*(24  hrs.) 

FLD - PF(K) - 


(B-17) 
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72  tt  *  recover cblo  dot*.  Biock  71  accumulate*  the  r>ert  ec  over at 

dose.  Block  72  teat*  for  aare  death.  Block  74  test*  the  day'*  do*e  vi.  )Uxla.r. 
ERD.  If  there  la  a  new  Maxirun  ERD,  it  1*  *et  in  block  75,  and  the  day  and  tine 
•re  incremented  in  blocks  76  and  77.  If  tKe  peak  ERD  (D3)  has  been  passed 
(D  <  D3),  then  the  next  point  of  interest  is  whether  or  not  the  allowable  level 
has  been  reached.  This  test  is  block  78.  If  the  dose  is  above  the  allowable, 
block  79  testa  to  aee  if  the  allowable  is  reachable--i.e. ,  nonrecoverable  dose 
is  less  than  the  allowable.  If  the  nonrecoverable  dose  exceeds  the  allowable, 
then  the  program  is  in  a  loop  as  the  length  of  stay  in  shelter  would  approach 
infinity — therefore,  the  error  message. 

Blocks  80  through  92  calculate  the  probability  of  fatality  and  the  probabili 
ty  of  noneffectiveness. 

The  preceding  explanation  should  give  some  understanding  of  the  computer 
program  and  how  it  represents  the  mathematical  formulation.  This  program,  then, 
generates  the  possible  data  combinations  that  are  analyzed  in  this  volume. 
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Code  Kumber  Series  1XX0 
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) 

SET 

Total  ERD 

D=0. 

Recoveral le  dose 

D1=0. 

Non-rec.  dose 

D2=0. 

Max.  ERD 

D3=0. 

Day  under 

calculation 

N=1 

Day  max.  ERD 

occurred 

Nl*l 

Day  ERD  return 

to  allowable 

N2=l 

Time  in  hrs. 

T=TA 

©r 


Increment  Time 
by  3  min. 

T  *  T  +  .  05 


ESH 

wmWss 

T3 

l _ 

_ 

©r 


RDN 

»  1.  -  COSF  (RDN)| 


RON 

=  R*RDN*T~ 

2.0 


*  D  +  RDN*.  05 
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Code  Number  Series  6XX0 
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PROGRAM  MAINLINE 

0 1  MENS  I  ON  PFN(IO)  „  _ _  _  .  __  . . 

N2-N2 

KODE-0  _  _ _ _  _ _ . . . . 

READ  INPUT  TAPE  5, lOod.RL ,RH,RS 

_ READ  INPUT  TAPE  5 .  1000 . TAL . TAH. TAS _ 

READ  INPUT  TAPE  5, 1000 , T2L , T2H, T2S 

REAOJNPUT  TAPE  5 , 1 OOOj  Q4L  i  D4H*  f)4J» _  _ 

READ  INPUT  TAPE  5 . 1000 , ZL . ZH , ZS 

READ _JNPUT  TAPE  5.  1000 .  FL  .  FH, FS _ _ _  _ 

READ  INPUT  TAPE  5.  1 000 . 8L  , BH,  BS 

_ READ  INPUT  TAPE  5.  1  OOP .  E TAL  .  E TAH.  ET AS _ ; _ _ 

READ  INPUT  TAPE  5 , 1001 , PFN ( 1 > , PFN I  2  1 1 PFN I  3  I , PEN  I  A ) , PFN ( 5 I , PFN I  6 ) , 

IPFNI 7),PPN(8l,PFNj9> _ _ 

READ  INPUT  TAPE  5, 1002,KMAX 

1000  FORMAT  (  3F 1 5  .5i ) _ 1 _ 1 _ 

1001  FORMAT ( 9F6.0) 

1002  FORMAT!  110) _ _ _ 

WRITE  OUTPUT  TAPE  6,8933 

8933  FORMAT!  115H _ R _ TA _ I2_ _ 04. _ Z  _ _ F  8 

1  ETA  K  PI  P2  D3  Nl  0  N  02  1 

2110  _R=RirWS  _ _ : _ J, . . . 

2120  IF(R-RH)2130, 2130, 8820 

2130  H-R+RS _ _ 

2140  TA  =  T  AL- TAS 

2150  IFITA-TAHI2l60j2160,2120_ _ _  _ 

2160  TA=  T  A  +  f  AS 

2170  _T2  =  T2L-T2S _ ... _ _ _ _ _ 

2180  IF  I T2-T2H)2190,2190,2150 

2190  T2  =  T2»T2S _ _ _ 

2200  D4-04L— D4S 

2210  IF (04-D4H) 2220, 2220, 2180 _ 

2220  D4-D4+04S 

2230  Z-ZL-ZS  _  _ _ _ _ 

2240  IF(Z-ZHI2250»2250»2210 

2250  Z*Z  +  Z  S _ _ 

F=FL-FS 

22  70_  IF  (F-FH) 2280, 2280, 2240 _ _ _ _ _ 

2280  F=F+fS 

'  ■ _ B-BL-BS _ j, _ _ _ 

2300  1F(B-BHJ2310, 2310, 2270 

2310  B-B  +  BS _ _ _ 

ETA=ETAL-FTAS 

2330  IF  !ETA-ETAH1 2340, 2340, 2300 _ _ 

2340  ET  A*  ETA  +  ETAS 

2350  DO  9000  K-l.KMAX _ _ _ _ _ __ 

3360  IFIZ-l. 1)3370, 3380, 3380 

3370  T3»( 1.1333+ETA1-TA _ _ 

GO  TO  4000 

_3380  IF !Z-1.  131 3390, 3400, 3400 _ _ 

3390  f  3  =  C 1 • 175+ETA 1 • TA 

GO  T0_4000  _ _  _ 

3400  T  3  * ( 1 • 2 6 ♦ E  T  A  )  *  T  A 

4000  CONTINUE _ 

4410  D-0  • 

pi-0.  _  _ _ _ _ _ 
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4420 

T-TA 

IF(T-(T2*TA) 144  30*  5480*5480 

4430 

T-U0.05 

4440 

RDN-(T-TA)«*3.i4l6/T3 

4450 

RON* l .O-COSF ( RDN ) 

_  4460 

RDN*(R«RDN)MT*»ZZ)/2.0 

4470 

0*0+ ( RDN-0.05 ) 

GO  TO  4420 

5480 

T-T+O.l 

5490 

RDN- (T-TA)»3«1416/T3 

5500 

RON-1. O-COSF(RDN) 

RDN*(R»RDN)'*(T-*ZZ)/2.0 

5520 

0*0* ( R0N*0.2/PFN( K ) ) 

5530 

1F(T-(TA->T3)  )5550. 5550. 5540 

5550 

T-T*0.2 

GO  TO  5490 

* 

5540 

T-T*0. 1 

IF ( T- ( T  A+24.0 ) J5570. 5610. 5610 

5570 

IF (Z- 1.0 15590* 5580, 5590 

5580 

ARG-24.0/T 

FLO-LOGF ( ARG ) 

GO  TO  5600 

5590 

ZZZ-l.O-Z 

FLO*( (T*»ZZZ)-< (TA*24.0)**Z ZZ))/(Z- 1.0) 

5600 

D-0+(FL0*R/PFN(K» ) 

5610 

IF (0-700.) 6000, 9000, 9000 

6000 

CONTINUE 

6640 

Dl-F-D 

6650 

02* ( 1 .0-F ) *0 

66  70 
6680 


6730 

6740 


6750 

6760 


6780 

6790 


7710 

7000 


7800 

TfllQ 


7820 

7830 


D3-0 

N»N»I 


6690 

6700 


6710 

HIQ 


02-02M  1  .pl-F  )  *FL0 

,D»gt*P2 


FLO* ( R*24  .'0 )«(T»»ZZ)/PFN(K) 
01»( l.-B)r01»(F»FLP) _ 


IF (0-700. 16740*9000*9000 
IF (D-03) 6780. 6750* 6750 
03-0 


Nl 


6770  T«T*24.0 


•  —  i  w 

GO  TO  6680 _ 

IF (D- 04) 7000 *6790, 6 790 
IF(D2-D4 >6770, 7710*7710 
KODE-1 


KODE-1 
CONTINUE 


N2-N 

1F(P3-1Q0«  >9000*7820*7820 


IF  1 03-200. I  7880*7830,7830 
IF (03-300, 17890, 7840, 7840 


7840  IF (03-700. ) 7900,9000,9000 

7880  P2  =  ( 03-  100  •  )  / 200  .  _  _  _ _ _  ______ 

7910  Pl=n. 

CO  TO  8000  _  _  . . . .  _ _ _ 

7890  P2=( ( 3.0*031-100. 1/1000. 

_ CO  TO  7920 _ _ _ 

7900  P2=C  700.-031/500. 

7920  Pl*(03-200.  1/500.  _  _ _ _  _ _  _ , 

8000  CONTINUE 

6950  MR  I  T  6  OUTPUT  TAPE  L , 6934j KOOE* R , T A , T 2 , 04 , Z t F » B , ET A , PFN ( K ) , 
1P1,P2,D3,N1,0,N2,D2 

6934  FORMATI I4.4X.F6.0.FB.1.F8.3.F6.0,F6.2.2X.F5.3.2X.F5.3.2X,F5.3, 
1F5.0,F7.3»F7.3»F7.1» 1 5 , F  7 , 1 ,  1 5, F9 .2  I 

9000  CONTINUE  _  _ _ _ _  _  . 

K0DE*0 

_  _G0  I0_2330  _ _ ' 

8820  CONTINUE 

" . ~ . END  SENSPROG  "  . .  " . — . — . — : . .  ; . ~  - . ; - - r . -- 


END 

FINIS 
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